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It is shown within the weak-coupling model that the macroscopic superconducting anisotropy for 
materials with the gap varying on the Fermi surface cannot be characterized by a single number, 
unlike the case of clean materials with isotropic gaps. For clean uniaxial materials, the anisotropy 
parameter 7(T) defined as the ratio of London penetration depths, X c /X a b, is evaluated for all T's. 
Within the two-gap model of MgB2, j{T) is an increasing function of T. 



A remarkable confirmation for the observed two-gap 
structuraJd of superconducting MgB2 came from solving 
the Eliashberg_muations for the gap distribution on the 
Fermi surface.QLj According to this, the gap on the four 
Fermi surface sheets of this material has two sharp max- 
ima: Ai ss 1.7 meV at the two 7r-bands and A2 ~ 7meV 
at the two er-bands. Within each of these groups, the 
spread of the gap values is small, and the gaps can be 
considered as constants, the ratio of which is nearly T 
independent. In this situation, a weak-coupling model 
with two gaps on two parts of the Fermi surface may 
prove useful in relating various macroscopic properties of 
MgB2. Starting with Ref. ||, the two-band models were 
studied by many, see, e.g., Ref. [9] and references therein. 
The focus of this work is on the macroscopic supercon- 
ducting anisotropy 7. To a large extent, motivation for 
this work was to understand why experiments on differ- 
ent samples of MgB2 done with-different techniques yield 
widely varying values for 7.E3il3 

The anisotropic Ginzburg-Landau (GL) equations, 
derived for clean superconductors with an arbitrary 
gap anisotropy in the-jseminal work by L.Gor'kov and 
T.Melik-BarkhudarovJlZI led to the commonly used con- 
cept of a single parameter 7 defined as £ a /£c = X c /X a (£ 
is the coherence length, A is the penetration depth, and 
a,c are principal crystal directions). Formally, this came 
out because the same mass tensor enters both the first 
GL equation which determines the anisotropy of £ (and 
of the upper critical fields H C 2) and the equation for the 
current which defines the anisotropy of A. However, it 
has been later shown by S.Pokrovsky and V.Pokrovsky 
in the work on the GL equations for anisotropic gaps in 
the presence of impurities, that 7 in fact depends on the 
impurity scattering, i.e., it might be sample dependcnt.ta 

In the following the near- T c result of Ref. |l| is repro- 
duced using the Eilenberger formalism. Moreover, the 
ratio X c /X a for arbitrary temperatures T is derived for 
the clean case. It is shown that for MgB2, X c /X a should 
increase with increasing T, the result which calls for ex- 
perimental verification. 

We begin with the quasiclassical version of— the BCS 
theory for a general anisotropic Fermi surface xB 



g 2 = 1 - // + , 
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A(r, v) = 2ttTN(0) £ (V(v, V )/(v' , r, u>)) v , . (4) 
j = -47r|e|iV(0)ri m ^(v 3 ). (5) 
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Here v is the Fermi velocity, IT = V + 27riA/0o; A is 
the gap function, /(r, v, w), / + , and g are Eilenberger 
Green's functions, N(0) is the total density of states at 
the Fermi level per one spin; fiu — ixT(2n + 1) with 
an integer n. Further, r is the scattering time on non- 
magnetic impurities and lud is the Debye frequency. The 
averages over the Fermi surface weighted with the local 
density of states oc l/|v| are defined as 



(X) = 



(27r) 3 ?iAf(0)|v| 



X . 
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Commonbi. the interaction V is assumed 
factorizable,EJ V(v,v') = V O(v)Q(v'), and one looks 
for A(r, T; v) = *(r, T) fi(v). Then, the self-consistency 
Eq. (Q) takes the form: 

um 

*(r, T) = 2ttTN(0)V £ (0(v)/(v, r, u)) . (7) 

The function Q (v) can be normalized by requiring 
that the critical temperature T c q for the clean material 
(r — > 00) is given by the standard isotropic weak-coupling 
model with the effective interaction VoEil 



(tt 2 > = 1 



(8) 



As usual, we incorporate T c o in the Eilenberger system 
using the identity 
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N(0)V " T c0 ^fhuj 
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Substitute this in Eq. (J?]) and replace u>r> with infinity 
due to the fast convergence: 



vll/ = 2Ag/Ti - 2luJ + (g(f) - f(g))/r , (1) 
-vIT*/ + - 2A*g/h - 2uf+ + (g(f+) - ! + { 9 ))/t, (2) 
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Effect of nonmagnetic impurities on T c . It is long 
known that scattering by nonmagnetic impurities, sup- 
press T c provided the gap is weakly anisotropic J23l23 The 
suppression is readily obtained from Eilenberger equa- 
tions without assuming that the anisotropy is weak. In 
zero field, all quantities are coordinate independent; be- 
sides, as T -> T c , g -► 1. Then, Eq. (Q) gives 



fiio' V 2lut ) hu>' 
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where u>' 
obtain: 



= u + l/2r. Substitute this in Eq. @ to 



w>0 
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Hence T c = T c o for r — > oo and any gap anisotropy; the 
same is true for the isotropic gap (f2 = 1) and any t. 
This equation can be written as 



ln^ 



(1 - (tt} 2 ) 



1 



/' 
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where /i = Ti/2t{T c t and tp is the di-gamma function. 
For a weak anisotropy (ft) 2 =j-L — \ with \ "C 1, this 
reduces to Hohenberg's resultE3 Formally, Eq. ( |13] ) is 
reminiscent of the case of magnetic impurities; the factor 
1 — (Q) 2 , however, makes a difference. For /i « 1, one 
has 



T r =T, 



cO 



8^ 
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For large fi's, unlike the case of the magnetic pair- 
breaking, we obtain 



T c = T c0 [A Q (0)T/h]( n y 



(15) 



where Ao(0) = 1.76 T c q. Therefore, T c does not turn 
zero at a finite r_nnless (ft) = as, e.g, for the d-wave 
superconductors.E£l 



Anisotropy near T c . As is seen from Eq. (|1 1Q , im- 
purities cause isotropization of amplitudes /, and one 
expects the macroscopic anisotropy to be suppressed by 
scattering. To address this question, one has to derive 
the GL equations iu the presence of impurities following 
basically the workt-3 for clean superconductors. As men- 
tioned above, the same mass tensor enters both the first 
and the second GL equations. We focus-.on the current 
equation because this is an easier task.El Within Eilen- 
berger formalism this is done in the clean case by expand- 
ing / near T c in two small parameters: A/Huj ~ ySt and 
vIIA/frw 2 - £oA/£T c ~ St (here St = (T c - T)/T c and 
£o is zero-T coherence length): 



^ Tiu) 
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Substituting this in Eq. ([!]) one obtains a = —1/2. The 
second GL equation follows by using Eq. ( p|) i n which we 
substitute jwl- //+/2 with /'s of Eq. <M): 
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In the London limit W = ^>oe l6 with a constant \&o, and 



.0«AV(V» + jjA),. 
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with 



(A 2 ): 
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The anisotropy parameter follows: 
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In fact, this is the result of Ref. [L7| 

In the presence of impurities, the first order term in 
the expansion ( |16 
defined in Eq. (Til 



should have the form (11). With D 
, we verify readily that 



D vUD 

~ hU ~ mJ 2 
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satisfies Eq. (hi). Writing D in the form: 



d = * in 



2lotJ 

we obtain with the help of Eqs. (JsJ) and 
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In the London limit, we have 
16n 2 e 2 N(0)T c 
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This is the result of Ref. |lj. In the clean limit it reduces 
to Eq. (p(i|), whereas the effect of impurities on 7 depends 
on the order parameter symmetry. 

For the d-wave symmetry (A) = and f2' = fi. In 
other words, the strong T c suppression notwithstand- 
ing, non-magnetic impurities do not affect 7. For order 
parameters with a non-zero (A), the strong scattering 
erases the effect of gap anisotropy on 7 altogether: 
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Hence, in the dirty limit, all parts of the Fermi surface 
contribute evenly to the anisotropy parameter as is the 
case for isotropic gaps. 

T dependence 0/7 = A c /A a - To address this question 
in the full temperature range one has to study weak su- 
percurrents, i.e., turn to Eq. (||). We consider only the 
clean case for which /o, go in the absence of currents are: 



fa — fo 



Ao 



9o 



B 1 



fe2 2 
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(27) 



in general, both Ao and /3 depend on . A weak super- 
current causes the order parameter A and the amplitudes 
/ to acquire an overall phase 6(r). We look for the per- 
turbed solutions in the form: 



A = A e ie , / = (/ + / 1 )e* 9 , 
/ + = (/o + /i + )e- <e , .9 = . 90 + .91, 



(28) 



where the subscript 1 denotes corrections. In the London 
limit, the only coordinate dependence is that-of the phase 
6, i.e., can be taken as r independent .cil Equations 

(Ul) then give: 



A 5i - hufi = ihf Q vP/2 , 
A o5 i - Juuf+ = ihfovP/2 , 
2.90.91 = -/o(/i + /i + ) ■ 



(29) 



Here, P = V6 + 27rA/0 o = 27r a/tftq. To evaluate the 
current (^J), one solves the system ( P^ ) for g\: 



.91 



A 2 



(30) 



Then one obtains the London relation between the cur- 
rent and the vector potential, 4irji/c = — (A 2 )^, 1 ^, with 



(A 2 ) 



//•• 



f3 3 



The anisotropy parameter now reads: 
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As T -► 0, we have 27TTA 2 , XL /?~ 3 1, and 



7 2 (0) 



(31) 



(32) 



(33) 



Note that the gap and its anisotropy do not enter this re- 
sult. The physical reason for this is in the Galilean invari- 
ance of the superfluid flow in the absence of scattering: 
all charged particles take part in the supercurrentEj 

Near T c , /?" 3 -> 7C(3)/8tt 3 T c 3 , and we obtain the 
GL result ( p0[ ) which amplifies contribution of the Fermi 
surface pieces with large gap to the parameter 7. Thus, 



the anisotropy parameter depends on T, the feature ab- 
sent in superconductors with isotropic gaps. 

It is worth noting that in literature the supercon- 
ducting anisotropy is commonly referred to as the ratio 
H C 2.a/ H C 2 iC , an important figure for applications, but a 
difficult quantity to evaluate for anisotropic Fermi sur- 
faces, not to speak about anisotropic gaps. On the other 
hand, measurements of the A-anisotropy such as, e.g., 
the high-field torque technique is more demanding than 
the resistive determination of H c2 since one has to work 
in the reversible domain. One should also note that, 
theoretically, the ratios of H c2 s and of A's are not nec- 
essarily the same, except near T c where their equality is 
provided by the GL theory. 

It is of interest to examine the consequences of oux-r^ 
suits for MgB 2 . The reported 7's vary from 1.7 to SjETt 3 ] 
or even higher as in Ref. [l6]. In all these reports, differ- 
ent techniques for extracting the anisotropy and samples 
with different resistivity ratios were used. 

Consider a model material with the gap anisotropy 
given by 



O = ni,2 , vef. 



1,2 : 



(34) 



where F\ , F 2 are two sheets of the Fermi surface. Denot- 
ing the densities of states on the two parts as Ai^, and 
assuming the quantity X being constant at each sheet, 
we obtain for the general averaging (|J): 

(X) = (XxNi + X 2 N 2 )/N{0) = v x X t + u 2 X 2 , (35) 

where we introduce normalized densities of state v\ t2 = 
Ai,2/A(0) for brevity. We have then instead of Eq. (@): 



ttlisi + Vl\v 2 = 1 . 



V\ + v 2 = 1 



(36) 



We also assume that the two parts of the Fermi surface 
have the symmetries of the total, e.g., (v)i =0 where 
the average is performed only over the first Fermi sheet. 
Within this model, Eq. (p2h reads: 



E^of 



* = 1,2, 



(37) 



where (i, = yVw 2 + V 2 (T)Of . 

Based on the band structure calculations, the relative 
densities of states V\ and v 2 of our model are « 0.56 and 
0.44.BB The ratio A 2 /Ai = O2/O1 « 4, if one takes 
the averages of 6.8 and 1.7 meV for the two groups of 
distributed gaps as calculated in Ref. [| Then, the nor- 
malization (f36) yields f^i = 0.36 and O2 = 1.45. 

Now, we have all parameters needed to solve the self- 
consistency equation ( [To| ) for ip(T) with / = A//3 (the 
clean case). This is done numerically and the result is 
shown in Fig. [l] along with A^(T). 
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FIG. 1. Temperature dependence of the two gaps calcu- 
lated as Ai — tp(T)Q,i versus T/T c . The upper curve is /S.2/T c , 
the lower one is Ai/T c , and the middle curve is ip(T)/T c eval- 
uated as described in the text. 

To evaluate j{T) of Eq. (^) we use the averages over 
separate Fermi sheets calculated in Rcf. = 33.2, 

- 42.2, («2) a = 23, and (v 2 c ) 2 = 0.5 x 10 14 cm 2 /s 2 
The result of numerical evaluation of "f(T) is shown in 
Fig. | 




FIG. 2. Temperature dependence of 7 = A c /A a versus T/T c 
for clean MgB2 calculated with parameters given in the text. 

The ratio of A's can be obtained, e.g., from the angular 
dependence of the reversible torque on single crystals in 
intermediate magnetic Mds tilted relative to the prin- 
cipal crystal directions nH Some torque data for MgB2 
were reported by M. Andst et alpj but the T depen- 
dence was not examined in detail. The ratio H c i,abl Uc2,c 
was shown to drop with increasing T from about 6 at 
15 K to 2.8 at 35 K (see also Refs. HIH). The authors 
estimate this ratio as w 2.3 — 2.7 at T c . Near T c , the ra- 
tios of H C 2S should coincide with the ratio of A's. In this 
work we estimate the latter as ~ 2.4, see Fig. §. Given 
this agreement and the prediction made here that the 
T dependence of 7 = X c /X a b should be opposite to that 
of the observed behavior of H C 2, a b/ H c2 , c , the detailed 
studies of 7 by torque or other methods are desirable. 
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